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Abstract 

We investigate the detection of a pseudoscalar cj) that couples electromagnetically via an 
interaction ^g(j)FF. In particular, we focus on the conversion of pseudoscalars produced in 
the sun’s interior in the presence of the sun’s external magnetic dipole held and sunspot- 
related magnetic helds. We hnd that the sunspot approach is superior. Measurements by 
the SXT on the Yohkoh satellite can measure the coupling constant down to g = 0.5-1 x 
10“^°GeV“^, provided the pseudoscalar mass m < 7x10“® eV, which makes it competitive 
with other astrophysical approaches. 
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I. INTRODUCTION 


Light neutral pseudoscalars arise naturally as a result of spontaneously broken global 
symmetries [^. They will be truly massless if the symmetry is not anomalous. Although 
some models can be directly tested in colliders, some of the best limits come from astro¬ 
physics. 

Such particles generically couple to the electromagnetic held through couplings of the 
form 


C = . ( 1 . 1 ) 

where F^^ is the dual of the electromagnetic held F^'^ and 5 ^ is a coupling constant. This 
coupling can lead to the interconversion of photons and pseudoscalars in a background 
electric or magnetic held. 

It is reasonable to attempt to use such a coupling to attempt to detect or exclude 
pseudoscalars of a given mass or coupling. There has been a recent resurgence of interest 
in detection using such couplings. Several astrophysical approaches have been considered 
for detecting these electromagnetic couplings 0 , including penetration of light from distant 
stars through intervening opaque clouds 0 , the ehect on evolution of stars and the 
appearance of X-rays generated in the cores of distant stars that have been double converted 
to pseudoscalars and back to photons again j^]. In all but the hrst of these methods, the 
pseudoscalars are generated from thermal photons in the presence of nuclear electric helds. 

The sun is the nearest and hence relatively brightest source of pseudoscalars. Naturally it 
has been the source of some recent pseudoscalars search proposals and experiments. Lazarus 
and his collaborators searched for solar pseudoscalars by attempting to detect a signal from 
pseudoscalar conversion in a static magnetic held on the earth [0. The probability of 
pseudoscalar conversion in such a helioscope is proportional to For this experiment, 

where the magnetic held is i? = 2 . 2 x 10 ^gauss and its length £ = 1 . 8 m, failure to detect a 
signal set the limit ai g = 3.6x10“® GeV“^, provided the pseudoscalar mass m < 0.03 eV. A 
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more expansive search is taking place in Novosibirsk, where a gimballed accelerator dipole 
magnet tracks the snn With B = 6xl0"^ganss and £ = 6m, the detector shonld be able 
to provide sensitivity down to abont g = 10“^°GeV“^, provided m < 0.1 eV. Given that 
maintaining a mnch stronger and larger magnetic held than the one used in Novosibirsk may 
not be possible, the limit from the Novosibirsk search may be the best attainable from this 
method of pseudoscalar-photon conversion on earth. 

There are, however, other sources of magnetic held in which solar pseudoscalars can 
convert to photons. On the surface of the sun, there exist many highly dynamic magnetic 
structures. Pseudoscalars emerging from the sun’s interior may convert to X-rays in the sun’s 
own external magnetic helds. These X-rays should be detectable by solar X-ray telescope, 
such as the SXT on the Yohkoh satellite. 

Two strong candidates for the source of pseudoscalar conversion in the sun’s exterior are 
the dipole held and sunspot-related helds. The sun’s dipole held, although rather weak, has 
an ehective length of a few hundred thousand kilometers. Sunspots have held strengths of 
a few kilogauss on the sun’s surface with an ehective length of no less than a few thousand 
kilometers. We hnd that the magnetic held of the sunspot produces a signihcantly larger 
X-ray count rate than the sun’s general dipole held. Furthermore, using the magnetic held 
of a large sunspot, it is possible with Yohkoh’s SXT to set the limit g < 0.5-1x10-1° GeV-i, 
provided m < 7xlO-°eV if no X-ray signature from pseudoscalar conversion is detected. 
This makes this method competitive with other approaches. However, the region we are able 
to exclude is almost completely covered by the horizontal branch star approach, which covers 
a higher mass range and couplings down to GxlO-n GeV”!. Because of the importance of 
redundant checking in astrophysics, and the relative ease of completing this analysis, we feel 
this approach is worthy of further consideration. 

In section 2 we discuss the production of pseudoscalars from the sun’s interior. In section 
3 we discuss their conversion to photons in the presence of the sun’s external magnetic helds 
and the limit which such conversion can achieve. Detection of converted photons by Yohkoh’s 
SXT is discussed in section 4. We summarize our hndings in section 5. 
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II. PRODUCTION OF PSEUDOSCALARS IN THE SUN 


Pseudoscalars with a coupling of the form (O) are produced primarily by the Primakoff 
process. In vacuum the cross section for this process is logarithmically inhnite, due to the 
infinite range of the Coulomb potential. However, at finite temperature and in a dense 
medium such as the sun’s, Debye screening cuts off this divergence and leads to a conversion 
rate per unit time for photon to pseudoscalar of P 


P = 


327r 


1 + 


K 

4 ^ 


In 1 + 


4^2 




( 2 . 1 ) 


where E is the energy of a photon, T is the temperature, and k is the inverse Debye screening 
length. Noting that all except heavy elements, which constitute less than 1% of the total 
mass fraction of sun’s interior, dissociate completely, we shall assume for simplicity the 
complete dissociation of all elements in the sun which gives 

dTTQ! 






( 2 . 2 ) 


where Hi and Zi are the number density and charge of each component of the local density. 


respectively, and a is the hne structure constant. Multiplying (|2.1|) times the density of 
thermal photons, we hnd the pseudoscalar production density per unit energy per second is 
given by 

dV 


(^2 ^ ^2^2) ^ - E 


(2.3) 


dVdE 87r3(e®/'^ - 1) 

where 

We now integrate this over a solar model, for which we use the model of Bahcall and 
Ulrich [^. The values of T and as functions of radius are graphed in Fig. 1. The nearness 
of the sun enables us to distinguish pseudoscalars emerging from different points on the sun’s 
disk. This suggests that we should not integrate pseudoscalar production throughout the 
volume of the sun, but only along a line of sight. Since the earth-sun distance is much 
greater than the sun’s radius, we shall treat the lines of sight as parallel. The total rate of 
production along a line of sight is thus given by 
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( 2 . 4 ) 


<ir g^-Rs, 

dx dy dE 

where the function F{E, h) is dehned as 


E{E,b), 


rVl — b'^ 

E(E,b) = / dz 
Jo 


(e^/T _ 1) L 


(E^ + ln(l + - E 


r=r?0\/z^+6^ 


, ( 2 . 5 ) 


where ^ and T are implicit functions of r, the distance from the center of the sun, and b is 
dehned so that bRQ is the impact parameter between the line of sight and the center of the 
sun. Values of E{E, b) for different energies at various b are graphed in Fig. 2. 

Note that most of the pseudoscalars are produced at 6 < 0.2. Because this is small, 
pseudoscalars emerge almost perpendicular to the sun’s surface. In subsequent calculations, 
we will treat the sun’s surface as hat, and the pseudoscalars as emerging directly along 
a perpendicular. This leads to errors proportional to 6^, or only a few percent. Once 
the pseudoscalars emerge from the surface of the sun, they can interact with the sun’s 
exterior magnetic held and back-convert into X-rays. We now proceed to hnd pseudoscalars’ 
contributions to solar X-rays. 


III. CONVERSION OF PSEUDOSCALARS IN THE SUN’S MAGNETIC FIELD 


In the presence of the sun’s exterior magnetic held, pseudoscalars can be converted into 
photons. The probability of conversion along a line of sight is given by O 


D(x,y) 


( 3 . 1 ) 


where we have introduced D{x,y) dehned as 

D{x,y)=[ B^{x,y,z)E^^^'> dz , (3.2) 

Jo 

where we have taken the line of sight to be the z direction, z = 0 being where the pseu¬ 
doscalars leave the sun, and z = L being the earth. The phase 9{z) has two sources: the 
mass of the pseudoscalar (if any) and the ehective plasma mass of the photon due to the 
presence of plasma outside the sun. In general 6{z) will be given by 
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(3.3) 


m = / 

JO 


^ (2'Jiane{z') m‘^\ ^ , 


h nieE 2E J 

where me is the mass of the electron and Ue is the total nnmber density of electrons, and m 
is the mass of the pseudoscalar. 

We will hnd later that there is only a signihcant effect if the pseudoscalar is very light. Let 
us assume for the moment that the pseudoscalar is massless or so light that its contribution 
to the phase is negligible. The rapidly falling nature of the sun’s coronal density outside the 
chromosphere assures that the phase will oscillate rapidly close to the sun and not change 
at all away from the sun. The rapid oscillations will wipe out the contribution from the 
magnetic helds very near the sun. For simplicity, we assume that the contribution is zero 
from this nearby zone, and that the phase 6{z) is constant outside this zone. We dehne the 
boundary height h by the somewhat arbitrary condition ^ = 9[oc) — 6{h). Using numerical 
models |]^ for the sun’s chromosphere and corona, we hnd h = 1200 km above the “surface” 
of the sun, dehned as the one optical depth level of the sun at a wavelength of 500 nm. We 
can replace the upper limit on integration by oo, because of the short range of the solar 
magnetic helds, so our formula for D is modihed to 


D{x,y)= B±{x,y,z)dz, 


(3.4) 


We then integrate the production rate (p.4|) and conversion probability (|3.1| ) over the 
regions of interest on the sun’s surface and divide by the solid angle of the sun to the earth 
to give 




dx / dy- 


dr 


-p = 




'0 


dxdyE{E,b) D{x,y) 


(3.6) 


dE AttL'‘J J ’dEdxdy 64ir*L2 
the diherential X-ray hux at earth due to solar pseudoscalar-photon coupling, where b is an 
implicit function of x and y. This general equation will be used to calculate the expected 
X-ray count rate coming from two components of sun’s exterior magnetic held; namely, the 
general dipole held, and sunspots. 

The sun’s general held is most apparent during the periods of sunspot minimum. During 
these periods, magnetographs measure a weak held of about one gauss around the poles 
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T^. Evidence that these polar helds have opposite polarities have led to the belief that 


this general held is a dipole. Although the solar dynamo model predicts the 22-year cycle 
of the solar dipole moment [|^, the exact behavior of the general held is still uncertain. 
For example, in 1957, just two years after the sunspot minimum, it was detected that the 
two poles had the same polarity for a period of 18 months. Nevertheless, we shall assume 
the existence of the solar dipole held and calculate its contribution to the X-ray count rate 
during periods of sunspot minimum. 

The height h where we can neglect plasma ehects is small compared to the scale of the 
dipole held, so to hnd D we begin integration at the sun’s surface. Since the magnetic held 
of a dipole at the equator is half the strength at the pole, and a dipole held falls as the 
cube of the distance, we can hnd the value of D along our line of sight to the center of 
the sun as \D(.\ = \BpRQ, where Bp is the magnetic held at the pole. We have explicitly 
calculated \D{x,y)\‘^ as a function of x and y, and we hnd there is very small variation over 
the region b < 0.2 where pseudoscalar production is important. So we can approximate 
\D{x,y) \ = \Dc\ everywhere, and hence obtain with the help of (^.5| ) 

d*hdip 




bF(E,b)dh . 


(3.6) 


dE 5127rV2 Jo 

Letting g = 10“^°GeV“^ and Bp = 1 gauss, values of for diherent energies are graphed 
in Fig. 3. As evident from Fig. 3, the sun’s dipole held as the source of the magnetic held 
for pseudoscalar conversion does not give a signihcant X-ray count rate. The signal would 
be detectable if there were no background, but this, sadly, is not the case. Fortunately, 
sunspots are a much more promising source, and we now consider them. 

As implied by (|0|) , strong magnetic held leads to a large pseudoscalar to photon conver¬ 
sion probability. On the sun’s surface, the strongest magnetic helds often arise from sunspots. 
Beside having strong helds, typically several kilogauss, other advantages of sunspots are that 
they are visible on the sun’s disk and that the X-ray background from their central dark 
region, the umbra, and their surrounding less dark region, the penumbra, is relatively small 
T^ . Although the shapes of sunspots are often highly irregular, it is convenient to model 
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sunspots as having a circular umbra surrounded by an annular penumbra. An approximate 


field distribution for this model is given by a gaussian |15], 


B^{r) = —rexp(-r ja ) , 


Tra^ 


(3.7) 


where is the component vertical to the solar surface, 0 is the total magnetic flux coming 
out of the sunspot, r is the radial distance from the sunspot center, and a is the characteristic 
size of the sunspot. The radius of the penumbra p is typically about = 2.1a^. Equation 
( ^T|) is the field roughly at z = 0. For the sunspots we are interested in, the size of the 
sunspot will be much greater than h = 1200 km where we can neglect plasma effects on the 
X-rays. Presumably, the fields change relatively little over this distance, and we will once 
again approximate h = 0. 

We now need to know the magnetic field, or, more to the point, D{x,y) due to the 
sunspot. Because the magnetic field is divergenceless, and there is no magnetic field at 
infinity, we can show 


V- D(x,y) = B^(x,y,z = 0) . 


(3.8) 


If we assume the magnetic field is rotationally symmetric above the surface of the sun, this 
is sufficient to relate D{r) to the magnetic flux 0(r) within a radius r, specifically 

(j){r)f 0[1 — exp(— 


D{r) = 


27rr 


27rr 


(3.9) 


It can be shown that if the assumption of rotational symmetry is dropped, this formula leads 
only to an underestimate of the effect. 

We now would like to calculate the total X-ray flux from a sunspot. The largest sunspots 
we will imagine will have a penumbra of size p = .O63i?0 |^. This is somewhat smaller than 
the scale on which our underlying pseudoscalar production changes. It therefore seems not 
unreasonable to approximate F{E, h) in ( |3.5| ) by the value at the center of the sunspot. Nu¬ 
merical calculations approximating the variation of E(E, h) throughout a very large sunspot 
show that this introduces at most a few percent error on the X-ray count rate. For a more 
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typical sunspot, p = .02 Rq, the error drops below 1%. With these two simplihcations and 
(B). we hnd that the differential X-ray ffux for sunspots is given by 


g^ReFjE^b) rp 
dE Jo 


27rr dr 


D{t] 


2 g^RQ(j)^F{E,b) 
1287r5L2 


[0.35] , 


(3.10) 


Letting = 10 ^°GeV ^ and optimistically setting 0 = 10^^ maxwells, which implies a very 
large sunspot with p ^ .06 Rq and i?z(0) ~ 3900 gauss, we find 


dA>s (F{E,h)' 

= 540 ' ^ ^ 


counts s ^cm ^keV 


dE 


keV5 


^-1 


(3.11) 


The maximum value at £' = 4keV and 6 = 0 gives about 410 counts cm“^s“^keV“^. 

Both (|375|) and (|3.10|) assume that pseudoscalars are massless particles. If they are 
not massless, the count rate given above can decrease significantly. As seen in (|3.3|) , the 
conversion of massive pseudoscalars along different points along the line of sight will produce 
photons of varying phase. For a sunspot with 0 = 10^^ maxwells, we find that pseudoscalars 
with mass m < 10“^eV and energy E = 4keV can cause a phase shift no greater than 7r/2; 
hence, the decrease in count rate would be minimal. However, for masses that are much 
greater, incoherent contributions of photons make the method proposed here unrealistic. 
Realistic pseudoscalar masses are either zero or typically in the range of 10“®-10“^eV [Q. 
For most of this range, the signal is reduced, probably to invisibility. But for masses near 
10“^ eV, there will be substantial signal, though somewhat reduced from the signal calculated 
in this paper. We will not calculate in the boundary region m ~ 10“^ eV explicitly, but will 
simply state that any limits obtained by this method using a detector that sees X-rays 
around 4keV will apply only for masses m < 10“^eV. 


IV. DETECTION OF X-RAYS CONVERTED FROM SOLAR PSEUDOSCALARS 

As evident from Fig. 2, most of the solar X-rays arising from pseudoscalar-photon con¬ 
version are in the energy range of 1-15 keV. At present, the best detector for finding such 
X-rays is the Yohkoh satellite, which is equipped with X-ray detecting instruments for the 
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primary purpose of studying solar flares. In particular, the Soft X-ray Telescope (SXT) on 
Yohkoh forms X-ray images in the 0.25 to 4.0 keV range 0- Because the expected signal 
detected by the SXT will be concentrated around 3 keV, we will only assume validity of our 
result for masses m < 7xlO“®eV. 

SXT has an effective area in this energy region of a little under 1 cm^. Multiplying SXT’s 
unhltered effective area with (|3.10|) and integrating through the energy range, we obtain the 
count rate detected by SXT from a sunspot due to pseudoscalars. Once again, assuming 
g = 10“^°GeV~^ and a very large sunspot such that 0 = 10^^ maxwells, we obtain SXT’s 
expected count rates for various b, as graphed in Fig. 4 with the maximum count rate being 
around 190 counts/s at 6 = 0. (In contrast, the signal from the whole quiet sun dipole held 
will be only 0.04 counts/s.) 

Unfortunately, the signal is dominated by background, which for such a large sunspot 
would be of the order of 4100 counts/s JT^. This does not necessarily make our signal 
unobservable. A single sunspot could be watched while it traverses across the face of the 
sun. The rate could be measured before it reaches the center of the sun, and compared to 
the rate while it is exactly on center. If the coupling g = 10“^°GeV“^, then there should be 
a four percent increase in the X-ray count rate. If systematic errors can be reduced below 
this level, then couplings at this level can be discovered or ruled out. 

If there are no systematic errors to worry about, then using a 10^ second exposure would 
allow one to see signals as low as 11 counts/s at the 5a level. Since the count rate scales as 
the fourth power of the coupling, this allows us to limit or discover the coupling down to 
g = 5x10“^^ GeV“^. This assumes we can reduce systematic problems to about the 10“^ 
level. 

Of course, there are systematic contributions. The penumbra may contain hot spots 
which produce many more X-rays than the rest of the sunspot. These hot spots may be 
transient, so that they might, for example, turn on coincidentally as the sunspot crosses the 
center. This effect can be limited by using the image to identify, isolate, and eliminate these 
hotspots. 
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Perhaps a more difficult problem is that the viewing angle of the sunspot changes as 

it crosses the face of the sun. This changes the apparent size of the sunspot, as well as 

changing the thickness of solar atmosphere through which it is viewed. One can argue that 
such effects should be proportional to 6^, the square of the impact parameter, and hence 
should be of the order of only a few percent or so as the sunspot moves from b = 0.2 to 

6 = 0. Thus it does not necessarily hamper us at the level of g = 10“^°GeV“^, and by 

correcting for it we can probably improve on this limit. 

There are a variety of other potential systematic headaches, such as angle-dependent 
detector efficiencies, evolution of sunspot as it moves, and perhaps others. None of these 
seem to be insurmountable, but they must all be systematically addressed before the value 
of g can be reliably discovered or ruled out. We believe that the ultimate limit from this 
approach using the Yohkoh detector will be in the range g = 0.5-1 x This is 

the limit for mass m < 7xlO“®eV. 


V. CONCLUSIONS 

We have found that searching for X-rays from the sun could limit the photon coupling of 
a light pseudoscalar. The X-ray production from the sun’s dipole moment is negligible, but 
with slightly optimistic assumptions about sunspots, we can obtain a limit using Yohkoh’s 
SXT from the umbral plus penumbral X-ray signal of about 

< 0.5-1 X 10“^° GeV“^ for m<7xlO“®eV. (5.1) 

This limit requires that we study the time evolution of the X-ray production from sunspots 
as they cross the center of the sun, looking for a brightening in the region b < 0.1. The 
exact limit depends to what extent systematic uncertainties can be eliminated. How does 
this compare with other approaches? 

Raffelt and Dearborn have previously placed limits on this coupling at the level 
10 “^'^ GeV“^ by considering the effects of pseudoscalars on the lifetime of horizontal branch 
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stars [^]. Naively, our limit looks as good or better than theirs, but only marginally so. 
In fact, their limit applies to masses as high as IkeV, and hence their limit applies in the 
all-important axion case, whereas ours does not. Furthermore, Raffelt has recently reex¬ 
amined the limit from this approach, and concluded that it can limit the coupling down to 
6 x10“^^ GeV“^ 1^, covering nearly all of the range of our detection. 

In a previous paper p, one of us considered what limit had been set or could be set 
by looking for X-rays from stars, particularly the nearby supergiant a-Ori. A limit of 
g < 2.5x10“^^ GeV“^ was set by reexamining old data, and it was pointed out that this 
limit could be improved to the level of 10“^^ GeV“^ if a dedicated observation were made. 
This is better than we believe the solar approach can achieve, but it applies only to masses 
below about 10-^°eV. 

Gombining the other two limits, we see that Yohkoh has a small window of opportunity 
for this method for masses in the range 10“^‘’-7xl0“®eV, if the coupling is very close to 
5x10“^^ GeV“^. How likely is this value? It is not favored by any model of which these 
authors are aware. For such a small coupling, eliminating systematic errors is crucial, and 
any “discovery” of a signal would be met with justihable skepticism. 

We therefore conclude that the approach we have proposed here is of only marginal value 
for investigating new portions of parameter space. On the other hand, more information 
is better. In the held of astrophysics, one never knows when some limit will be called into 
question, and it is always good to have redundant approaches to particle physics limits. 
Redundant approaches are as close as we can come to repeatability of experiments in astro- 
particle physics, so we should seize every opportunity we can hnd. Who knows? Maybe this 
approach will hnd evidence for pseudoscalars, or perhaps just lead to a better understanding 
of X-rays from sunspots. 
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Figure Captions 


Fig. 1. The temperature (solid line) and (dashed line) as functions of radius r based on 
the model of Bahcall and Ulrich [jTO . 


Fig. 2. The function F{E, b) for b from b = 0.0 (top curve) to 6 = 0.3 (bottom curve). The 
solid lines are at intervals of Ab = 0.1, and the dashed lines are at intervals of Ab = 0.02. 


Fig. 3. The expected X-ray signal due to pseudoscalar-photon conversion in the sun’s 
general dipole field, taking g = 10“^°GeV“^ and Bp = 1 gauss. 

Fig. 4. The expected X-ray count rate detected by Yohkoh’s unhltered SXT due 
to pseudoscalar-photon conversion in the magnetic held of a sunspot, taking g = 
10“^*^ GeV“^ and 0 = 10^^ maxwells. The sunspot’s center dehnes the impact parameter 
bR^. 
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